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The thawing quintessence model with a nearly flat potential provides a natural mechanism to 
produce an equation of state parameter, w, close to —1 today. We examine the behavior of such 
models for the case in which the potential satisfies the slow roll conditions: \{l/V){dV / dfj))^ <^ 1 
and {l/V){d'^V/d(f>^) <^ 1, and we derive the analog of the slow-roll approximation for the case in 
which both matter and a scalar field contribute to the density. We show that in this limit, all such 
models converge to a unique relation between 1 + vu, fl^, and the initial value of {l/V){dV/d(j>). We 
derive this relation, and use it to determine the corresponding expression for w{a), which depends 
only on the present-day values for w and Q^. For a variety of potentials, our limiting expression for 
10(0) is typically accurate to within Sw ^ 0.005 for w < —0.9. For redshift z ^ 1, Kj(a) is well-fit by 
the Chevallier-Polarski-Linder parametrization, in which w{a) is a linear function of a. 



I. INTRODUCTION 

Observational evidence [J, [3| indicates that roughly 
70% of the energy density in the universe is in the form 
of an exotic, negative-pressure component, dubbed dark 
energy. (See Ref. Q for a recent review). The observa- 
tional bounds on the properties of the dark energy have 
continued to tighten. Taking w to be the ratio of pressure 
to density for the dark energy: 



w = Pde/p 



DE, 



(1) 



recent observational constraints are typically 
— 1.1 !S w ^ — 0.9 when w is assumed constant 
(see, e.g., [1, Hf and references therein). 

We can consider two possibilities. If the measured 
value of w continues to converge to a value arbitrarily 
close to —1, then it is most reasonable to assume a cos- 
mological constant (a conclusion supported by both the 
Akaike information criterion 6] and common sense) . On 
the other hand, it is conceivable that the observations 
will converge on a value of w very close to, but not ex- 
actly equal to, —1. In this case, we must consider how 
such a dark energy equation of state might arise. 

One possibility, dubbed quintessence, is a model in 
which the dark energy arises from a scalar field [7, 
i, i, [HI- Caldwell and Linder [H showed that 
quintessence models in which the scalar field potential 
asymptotically approaches zero can be divided naturally 
into two categories, which they dubbed "freezing" and 
"thawing" models, with quite different behaviors. Thaw- 
ing models have a value of w which begins near —1 and 
increases with time, while freezing models have a value of 
w which decreases with time, with an asymptotic value 
that depends on the shape of the potential. (If the ob- 
servations converge to a value of w less than —1, more 
exotic models must be considered. We will not consider 
this possibility here). 

Thawing models with a nearly flat potential provide 
a natural way to produce a value of w that is close to. 



but not exactly equal to —1, since the field begins with 
w ~ — 1, and w increases only slightly up to the present. 
Furthermore, with a nearly flat potential and w ~ ~1, 
the dynamics of quintessence are considerably simplified. 
In this paper, we show that all such models converge 
to a single unique evolution. In addition to providing a 
plausible model for w near —1, such models (since they 
all converge to a single type of behavior) can serve as 
a useful set of fiducial models that can be compared to 
ACDM. 

In the next section, we reexamine thawing quintessence 
models and outline the arguments for these models. In 
Sec. 3, we present a general description of the behavior of 
such models. Our main results are given in equations (j23p 
and In Sec. 4, we discuss the arguments against 

thawing quintessence models. Our conclusions are sum- 
marized briefiy in Sec. 5. 



II. THE CASE FOR THAWING 
QUINTESSENCE 

We will assume that the dark energy is provided by a 
minimally-coupled scalar field, 0, with equation of mo- 
tion given by 



where the Hubble parameter H is given by 



H 



(2) 



(3) 



Here a is the scale factor, p is the total density, and we 
work in units for which SttG = 1. Equation ^ indicates 
that the field rolls downhill in the potential V{(j)), but its 
motion is damped by a term proportional to H. 

The pressure and density of the scalar field are given 



2 



by 



and 



P - 



(4) 



(5) 



respectively, and the equation of state parameter, w, is 
given by equation ([1]). 

As noted above, observations suggest a value of w near 
— 1. Caldwell and Linder [T^] noted there are basically 
two ways to achieve such a result (see also Ref. 
for a more detailed discussion of some of these issues). 
In the first case, freezing potentials, the field is initially 
rolling down the potential with w 7^ — 1, but it slows 
with time, driving w toward —1. (As emphasized in Ref. 
[1^, the tracking models introduced in Ref. are a 
subset of freezing models, but not all freezing models 
display tracking behavior). In thawing models, the field 
is initially nearly frozen at some value 4>q, with w = —1. 
Then as H decreases, the field rolls down the potential, 
and w increases with time. 

It is possible to produce w near — 1 at the present using 
either freezing or thawing models. Chongchitnan and 
Efstathiou |14| used Monte Carlo simulations to derive 
a set of potentials that yield w very close to —1 today. 
They found two classes of acceptable solutions: very flat 
potentials, and models in which the field evolves from a 
region with very steep slope in the potential to a region 
in which the potential is roughly flat. While neither type 
of model can be ruled out, we feel that the models with 
a nearly flat potential are clearly a more natural way to 
produce the desired present-day value of w near —1. 

A similar naturalness issue was raised by Bludman , 
who argued that the only way to achieve tracking mod- 
els with w near —1 today was for the model to contain 
a sharp change in the curvature of the potential at the 
present. Thus, we are presented with a double coinci- 
dence problem: why should the field be entering this spe- 
cial region of the potential at the same time that the dark 
energy is coming to dominate the matter, and why are 
both of these happening right now? In thawing models 
with a nearly flat potential, on the other hand, w never 
deviates very far from —1. 

Another argument in favor of these thawing models is 
that we already have strong evidence that the universe 
at one time underwent a period of vacuum energy domi- 
nation (inflation). Many models for inflation correspond 
to the sort of thawing models examined here the 
scalar field initially has w = — 1, but then rolls downhill 
to terminate inflation. 



Finally, Griest [17[ has suggested a solution to the co- 
incidence problem involving thawing fields. In his model, 
the universe contains a variety of scalar fields with vari- 
ous initial energy scales T/i(0o): * = 1; 2, 3, .... As the mat- 
ter or radiation density drops below a given ^^(^o), the 



universe undergoes a period of dark energy domination, 
but the field then thaws and slides down the potential, 
allowing matter or radiation to dominate again. Given 
enough of these fields, it would not be suprising to find 
ourselves in an epoch in which one of them is just be- 
ginning to dominate at present (see also the somewhat 
different model of Ref. p^V 

None of these arguments proves, of course, that a uni- 
verse with w close to —1 (but not equal to —1) must 
involve a thawing quintessence field with a nearly flat 
potential. However, they do indicate that such models 
are worthy of further study. 



III. EVOLUTION OF THAWING 
QUINTESSENCE WITH NEARLY FLAT 
POTENTIAL 

We will assume a scalar field with initial value <j}Q in 
a nearly flat potential V{(j)). Speciflcally, we will assume 
that at = 00, the fleld satisfles the slow-roll conditions: 



1 dV\ 



and 



1^ 

V df' 



< 1. 



(6) 



(7) 



The latter condition corresponds to a mass scale <^ 
^y((/)o). Setting V{4}q) roughly equal to the dark energy 
density at the present, we get <C 10"^'^ eV. This is the 
same (unnaturally small) mass that occurs generically in 
quintessence models for dark energy. 

In analyzing models for inflation, it is usually assumed 
that the scalar field dominates the expansion, and that 
V{(j)) » 0^/2. With these assumptions, along with the 
flatness of the potential given by equations ^ and ([7]), 
it can be shown that the term in equation ^ can be 
neglected, yiel ding the simple equation 3H(j) — —dV/dcj) 
(see, e.g., Ref. [ia|)- This is called the slow-roll approxi- 
mation. 

It is well-known that the slow-roll approximation fails 
for the case of quintessence [H, [H, • The basic rea- 
son is that the slow-roll approximation requires the scalar 
field to dominate the expansion. However, this is never 
the case for quintessence, since matter always contributes 
significantly to the total density. However, nothing pre- 
vents us from assuming the slow roll conditions on the 
potential (equations [S] and [7]) along with the requirement 
that w be close to —I today. Effectively, we are deriv- 
ing the analog of the slow-roll approximation for the case 
where the expansion is not dominated by the scalar fleld. 

At the late times which are of interest here, the uni- 
verse is dominated by dark energy (assumed to arise 
from a scalar field) and nonrelativistic matter; we can 
neglect the radiation component. We assume a fiat uni- 
verse containing only matter and a scalar field, so that 
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CI^ + Qm = 1- Then equations ^ and Q can be rewrit- 
ten in terms of the variables x, y, and A, defined by 



X = (b'/V6, 

y = VV{^)/3H^, 
1 dV 



A 



(8) 
(9) 

(10) 



and the prime wiU always denote the derivative with 
respect to In a: e.g., (f)' = a{d4)/da). (This discussion, 
from equation [8] through equation [TBI is taken from Refs. 

Then x^ gives the contribution of the kinetic energy of 
the scalar field to fi^, and gives the contribution of 
the potential energy, so that 



while the equation of state is 

2x^ 

7 = 1 + w = 



x'^ + y-^ 



(11) 



(12) 



It is convenient to work in terms of 7, since we are in- 
terested in models for which w is near —1, so 7 is near 
zero, and we can then expand quantities of interest to 
lowest order in 7. Equations ([2]) and Q, in a universe 
containing only matter and a scalar field, become 

x' = -^3x + X^y^ + ^x[l + x^-y% (13) 



where 



y' = ->^\J^xy+^y[l + x'^ -y'^], 
A' = -V6\^{T-l)x, 



d(j)^ ' \ d<p) 



(14) 
(15) 

(16) 



We now rewrite these equations, changing the dependent 
variables from x and y to the observable quantities U,^ 
and 7 given by equations pTjl and To make this 

transformation, we assume that x' > 0; our results gen- 
eralize trivially to the opposite case. We obtain: 



7' 
0.' 



-37(2-7) + A(2-7)v/3^, 
3(l-7)r!^(l-r!^), 

-\/3A2(r- 1)7^. 



(17) 
(18) 

(19) 



Note that equation (fT7|) also follows, in a trivial way, 
from the expression for w' given in Ref . [l^ . Finally, we 
will see that the equations simplify if we transform our 
dependent variable from a to Q,^{a). This gives us 



(^7 

dU,(k 



-37(2-7) + A(2-7)73^ 
3(l-7)r!^(l-ri0) 



(20) 



This change of variables is valid only if i}^ is a monotonic 
function of the scale factor; it breaks down at any point 
where dflc/j/da — 0. This condition is satisfied for most 
quintessence models and for all of the models we consider 
here; it is not satisfied, for example, in models in which 
oscillates in time |18l |. 
Equations (fT9|) and (f20|) are an exact description of 
the scalar field evolution for a:' > 0, but they do not 
yield any simple solution. At this point, we make two 
assumptions. Our first assumption is that 7 <|; 1, corre- 
sponding to w near —1. The second assumption is that 
A is approximately constant, so that 



X^Xo^-{l/V)idV/d(b) 



(21) 



= 00 



i.e., Aq is the value of A at the initial value of the scalar 
field 4>o before it begins to roll down the potential. Equa- 
tion ((2T|) follows from the slow-roll conditions, equations 
([6]) and ([7]), as we will show at the end of this calculation. 
Replacing A with Aq and retaining terms to lowest order 
in 7 in equation ([20|) yields the following: 



d'y 

dflrh 




(22) 



This equation can be transformed into a linear differential 
equation with the change of variables = 7, and the 
resulting equation can be solved exactly. For the models 
of interest here, we have the boundary condition 7 = 
at rt^ = 0. The resulting solution (reexpressed in terms 
of w) is 




Equation (j23p . along with the corresponding result for 
w{a) derived below, is our main result. It shows that for 
sufficiently flat potentials, all thawing quintessence mod- 
els with w near —1 approach a single generic behavior, 
with w{a) determined entirely by ^^(a) and the (con- 
stant) initial value of {l/V)(dV/d(j)). A graph of this 
generic relationship between w, Q^, and Aq is given in 
Fig. 1. 

Equation ([23]) shows that 1 -I- w - O(A^). Thus, 
our first slow-roll condition (equation [6|) insures that 
1 + w <C 1, as desired. The condition that A be nearly 
constant up to the present day can be quantified by re- 
quiring I A'/ A I <C 1. Taking 7 to be of order A^ in equation 
P^ . we obtain the condition 



1 d^V 
V~d^ 



1 dvy 



(24) 



The two slow-roll conditions, taken together, insure that 
this condition is satisfied. 
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FIG. 1: The value of (1 + TO)/Ag as a function of O,^ in thawing 
quintessence models with a nearly flat potential and w near 
~1. Here w and Q,^ are functions of the scale factor a, while 
Ao is the (constant) initial value of —{l/V){dV/d(l>). 
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FIG. 2: A comparison between ■w{Q,p) for V{(t>) = 0^ (dotted 
curve), V{(f)) = (j)-'^ (dashed curve) and our analytic result 
for w{fl4,) (solid curve) for (top to bottom) Ao = 1, Ao = 2/3, 
and Ao = 1/2. 



A sufficiently accurate determination of the present- 
day values of w and fi^ {wq and fl^o, respectively) 
uniquely determines the value of Aq for these models. 
For example, for fl^Q = 0.7 and wq = —0.9, we obtain 
Ao = 0.8. 

The way in which arbitrary potentials satisfying the 
slow-roll conditions converge to equation P5|) is illus- 
trated in Fig. 2. In this figure, the solid curve gives 
the behavior for ^(ri^) predicted by equation ((23)) . while 
the dotted and dashed curves give the true evolution for 
the potentials V — (j)^ and V = respectively, where 
we choose the initial value of (p such that Aq = 1, 2/3, 
1/2. As expected, agreement is poor for Aq = 1 and im- 
proves for smaller values of Aq. Since we dropped terms 
of order I + w in deriving equation (|23p , we expect the 
fractional error in 1 -I- w to be on the order of 1 -I- w. For 
w < —0.9, this translates into an error in w of dw ^ 0.01, 
which is apparent in Fig. 2. 

The behavior of the potential demonstrates an im- 
portant point: while negative power law potentials usu- 
ally give rise to "freezing" models 0, [13) El > they can be 
made to act as thawing models by an appropriate choice 
of (pQ. For example, when V{(p) — (p^", if (po ^ n, then 
equations ^ and ([7]) are satisfied, and the model behaves 
like a thawing model. This shows that any potential can 
give rise to the type of models discussed here, as long as 
V{(j)) has a region over which equations ([6]) and ([7]) apply. 

We can use equation (|18p to solve for fl^ as a function 



of a and thus determine w{a). Taking the limit 7 ^ 1 in 
equation (fT8|) gives fl'^ 



3flcj,{l — r2</,), with solution 



1 



"00 



1 a 



-3 



(25) 



where il^o is the present-day value of Q^, and we take 
a = 1 at the present. Equation (P5|) is identical to the 
expression for J^a as a function of a in the ACDM model, 
which is not surprising, as we are taking w near —1 (see 
also Ref. [1^). Equations (pS)) and (^5)1 together give an 
explicit expression for w as a function of a. Assuming 
a particular value wq for the present-day value of w, we 
can then eliminate Aq from this expression, so that w{a) 
is a unique function of wq and fi^o- We obtain: 



l + w = (l + wo) 
-(0,J - l)a"'''tanh"^ 




(26) 



A graph of w{a) is shown in Fig. 3 for several values of wq 
and rJ^o- Note that w{a) depends primarily on wq and is 
not very sensitive to the value of H^q for 0.6 < il^o < 0.8. 
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FIG. 3: Our analytic result for the evolution of w as a function 
of the scale factor, a, normalized to o = 1 at the present, in 
thawing quintessence models with a nearly flat potential and 
w near —1, for fl^o = 0.8 (dashed), fl^o — 0.7 (solid) and 
i^4>o = 0.6 (dotted). Upper three curves are for wo — —0.9; 
lower three curves are for wo ~ —0.95. 



Also, w{a) is a nearly linear function of a for a between 
0.5 and 1 (redshift z ^ 1), although the linear behavior 
breaks down for z > 1. The z ^ 1 behavior agrees well 
with the Chevallier-Polarski-Linder parametrization (2^ . 
f25j . in which w{a) is taken to have the form w{a) = wq + 
Wq(1 — a), although equation (|26p does not provide any 
particular insight into the origin of this linear behavior. 
However, in our case Wa is not a free parameter, but 
depends on wq and fl,po; for a fixed value of fJ^o, equation 
(|26p corresponds to a one-parameter family of models. 
For instance, for ri0o = 0-7, the linear fit to equation 
(PS)) for z £ 1 is roughly w — wq — 1.5(1 + wo){l — a), so 
that Wa ~ —1.5(1 + Wo)- 

Now we can compare the generic behavior predicted by 
equation (|26p with the actual scalar field evolution. In 
Fig. 4, we show this predicted behavior, along with w{a) 
for the potentials V — (f>'^,V — and V = exp(— A0). 
The value of (j>o for the power law potentials is chosen 
to give fl^Q — 0.7 and wq = —0.9. For the exponential 
potential, O^g = 0.7 and wq = —0.9 are fixed by the value 
of A. The typical errors here are Sw ^ 0.005, showing 
strong agreement between the true evolution of w{a) and 
our analytic expression for w{a). The error decreases as 
Wo decreases, so Fig. 4 gives an upper limit on the error 
in our approximation for wq < —0.9. 

In Fig. 5, we compare our limiting behavior for w{a) in 




a 



FIG. 4: The evolution of w as a function of the scale fac- 
tor, a, normalized to a = 1 at the present, with = 0.7 
and Wo — —0.9. Solid curve is our analytic result for the 
behavior of thawing models with a nearly flat potential and 
w near —1. Other curves give the true evolution for the po- 
tentials V{4>) — 4>^ (dotted), V((j)) = (^~^ (short dash), and 
V{(j>) = exp(— A0) (long dash). Dot-dash curve is a model 
with V{<l>) = exp(— A()i) in which A changes discontinuously. 

equation l|26|) to the SNIa observations. The likelihoods 
were constructed using the 60 Essence supernovae, 57 
SNLS (Supernova Legacy Survey) and 45 nearby super- 
novae, and the new data release of 30 SNe la detected 
by HST and classified as the Gold sample by Riess et al. 
[2, HEl ■ The combined dataset can be found in Ref. Q . 
It is clear that current observations do not exclude the 
thawing quintessence models we have considered here, 
although the observations are also obviously consistent 
with a cosmological constant. The maximum likelihood 
point actually lies below wq — —1, but we have not ex- 
tended our graph to wq < —1, as our derivation of equa- 
tion (|26p assumes a standard quintessence model with 
w > —1 at all times. 

We have shown that the slow-roll conditions, equations 
(O and ([7]), are sufficient to allow a thawing model to 
produce w near —1 today, but are they also necessary 
conditions? The answer is no, although violating these 
bounds with a thawing model requires rather unusual 
potentials. Following Ref. [l3|, we have produced a po- 
tential by joining the functions V[cf)) = exp(— Ai0) and 
V[(j)) = exp(— A2^) at (/) = 0. The potential is then con- 
tinuous, but A varies discontinously, so equation ([7]) is 
violated. In Fig. 4, we show w{a) for this thawing po- 
tential (dot-dash curve), where we have taken Ai = 0.2, 
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FIG. 5: The 1 — a (solid) and 2 — a (dashed) contours in the 
plane defined by the present-day values of il^ and w for the 
quintessence field (denoted fi^o and wq respectively), for the 
thawing behavior given by equation (|26|) . 



A2 = 1.04, and we have chosen ^0 to give wo — —0.9 
and Q^Q = 0.7. The evolution of (f> is clearly always in 
the thawing regime (dw/da > 0), and it gives w near 
— 1, but it docs not produce a functional form for w{a) 
resembling that of a nearly flat potential. Of course, this 
form for the potential is rather pathological (note that it 
also violates the upper bound on w' as a function of w 
for thawing models postulated in Ref. [13 )• 

One of the best-motivated thawing quintessence mod- 
els is the Pseudo-Nambu Goldstone Boson (PNGB) 
model [13]. (For a recent discussion, see Ref. [ll] and 
references therein). This model is characterized by the 
potential 



F(0)=M4[cos(</)//)-f 1], 



(27) 



and 00 can be taken to lie between and nf. Then the 
evolution of this model is a function of M, f, and ho- 
using equations (O and ([7]), we see that the slow-roll 
conditions are satisfied for all (j)o ii f > 1, while they 
cannot be satisfied for any 0o if / < 1. The latter result 
follows from the trigonometic function in the PNGB po- 
tential: its second derivative is large whenever the first 
derivative is small, and vice- versa. None of these results 
depend on the value of M. Thus, our results are a good 
approximation to the behavior of the PNGB model for 
the case where / > 1. 

Now we consider some related approximation schemes. 
Crittenden et al. [2^ analyzed quintessence models 
with w near —1 in terms of the parameter «(<?!)), defined 
through the equation 



dV/d(j) 



v{i + (t>/iH4)) 



(28) 



With w near —1, they took the evolution for to be 
given by equation (P5)) . and they approximated the evo- 
lution of linear function of 6: 



k{4>) = Kq + Kl{4> - 4>q). 



(29) 



Thus, the model of Ref. [23| has two free parameters, 
kq and Ki, which determine w{a). It is straightfoward to 
derive the equivalent of our equations ([23|) and ([26|) . For 
w as a function of ri^,, we obtain: 



1 



1 - n 



<i>Q 



4ki/3 



(30) 



Our corresponding result (equation [53]) is clearly distinct 
from this result for all values of kq and ki. Note further 
that equation (|30p implies 1 + wq — {2/3)KQil^Q, so we 
can express w as a function of wq, ^l,po, and ki alone 
(corresponding to equation [55)1 : kq drops out of this ex- 
pression: 



1 + w = (l + wo)a^ [n^QO^ + (1 - ^^0)] 



-(4ki+3)/3 



(31) 



Again, this result is distinct from equation (I26p , although 
the two expressions obviously can be made to converge to 
similar behavior by the appropriate choice of ki, since ki 
can be chosen to give good agreement with the exact evo- 
lution [2^. The main difference between this approach 
and ours is that our final result for w{a) contains no free 
parameters; it is a function only of wq and f2<^0j while the 
expression for w(a) in the form of equation ()31|) contains 
the fitting parameter ki. 

Neupane and Scherer [l^ considered the consequences 
of fixing X (as defined in equation[8]) to be a constant, a. 
With X = a, their relation corresponding to our equation 
(1231) is 



1 



a 

STL 



(32) 



While it is possible in such models to produce w close to 
— 1 at the present, equation ((32)) shows that these w « — 1 
models always act as freezing models, since increasing fl^ 
corresponds to decreasing 1 -\- w. 

One might argue that the correct "generic" model for a 
nearly flat potential should be a linear potential, i.e., con- 
stant dV/d(j), rather than constant {l/V){dV/d(j)). Linear 
potentials have been investigated previously by a num- 
ber of authors [Hil [M [H, [H, [Hi . Exact solutions 
for the linear potential have been derived for the case 
where the universe is scalar field dominated [s^l or when 
it is matter-dominated [Toj , but not for the intermedi- 
ate case. However, it is possible to use the techniques 
discussed here to derive an approximate solution. If we 
take 



V = Vo~ a^, 
then equation ^ has the solution 



/■*/ 


\ a{t) 1 


lt=ti 


Uitf)\ 



1 3 



dt. 



(33) 



(34) 
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The integrand can be reexpressed in terms of a and H{a) 
to give 



a 



a=a, H{a) \af 



da 



(35) 



This sohition is as yet exact. Now we make essentially the 
same approximations that we used earlier for 1 + w <C 1. 
We take H{a) to be given by equation but in deter- 
mining the total value of p, we approximate as a con- 
stant, and we take 1+w to be given by 1+w « (/)^/F((/)o). 
With these approximations, equation ([55)1 can be used to 
derive ^(ri^). The resulting expression for w{Q^) is iden- 
tical to equation (j23p . Furthermore, numerical integra- 
tion for the linear potential gives results in good agree- 
ment with our (wo^Wa) fit discussed above (see Fig. 4 of 
Ref. [31]). This supports the conclusion that our results 
feguations 1231 and I26 p represent a generic asymptotic be- 
havior. Of course, we cannot rule out the possibility of a 
more exact solution for the linear potential than the one 
we have outlined here. 



IV. THE CASE AGAINST THAWING 
QUINTESSENCE 

Now consider the arguments against the models con- 
sidered here. To avoid unknown quantum gravity effects, 
it is desirable for the energy scale of the scalar field to be 
below the Planck mass (unity in our units). Requiring 
(f) < 1 does not constrain the models presented here, as 
the potential can be modified to shift the value of (p to 
any desired value. In Ref. it was suggested that 

a possible constraint is < 1- Obviously, if this 

constraint is enforced, then all of the models considered 
here are ruled out, since we have only considered models 
with V' /V < 1 (equation [S]) . The implications of this 
proposed constraint are explored further in Ref. [l^ . In 
Ref. [1^, it was argued that the correct constraint is 
actually At/) < 1, where Acj) is the change in the value 
of (j) from its initial value to the present. Our models do 
satisfy this second constraint. 

Linder [l^ has noted that thawing models with Aq ^ 
1 occupy only a very small fraction of the phase plane 
defined by w and w' . This is certainly true; if one requires 
w to be very close to —1 at present, and assigns equal a 
priori weight to all phase trajectories in the w — w' plane, 
then models with Aq <C 1 are very unlikely. 

Huterer and Peiris (36j performed Monte Carlo simu- 
lations of quintessence models, sampling low-order poly- 
nomial potentials. (See also the related work in Refs. 
[13, They found essentially no acceptable thawing 

models were generated using this procedure. This is not 
surprising, since their procedure samples a uniform dis- 
tribution in the initial values of both [{l/V){dV/d(p)]'^ 
and {l/V)(d'^V/d(j)^), while our models require the ini- 



tial values of both of these quantities to be much smaller 
than unity. 

We do not dispute the conclusions in cither Refs. [l^ 
or [36 | ; they simply represent a different approach to de- 
termining the most plausible models. The models pre- 
sented here require a very flat potential. One can argue 
that the special nature of the potential makes such mod- 
els unlikely; however, we believe that an observed value 
of w near — 1 argues in favor of choosing such special po- 
tentials, while the fact that all such models converge to 
a similar evolution makes these models more interesting. 
These thawing models do require a fine-tuning of ^o; it 
must be chosen so that T^(0o) is approximately equal to 
the dark energy density today. 

The most serious problem with the models considered 
here is that there is currently no compelling observational 
evidence to favor them over a cosmological constant, as 
Fig. 5 shows. On the other hand, current observations 
do not rule out these thawing models. 



V. CONCLUSIONS 

Thawing models with potentials that satisfy the slow- 
roll conditions provide a natural way to produce w near 
— 1, and they all converge to a single, universal behav- 
ior. Such models are, in some ways, the opposite of the 
tracker models proposed in Ref. flrl . The tracker models 
are insensitive to the initial conditions, but they depend 
sensitively on the shape of the potential over the entire 
range of evolution of (j). The models discussed here, in 
contrast, depend only on the initial conditions, i.e., the 
value of V and its derivatives at 0o, but are insensitive 
to the shape of the rest of the potential. This situation 
arises because the field never rolls very far along the po- 
tential, and so never has a chance to "see" the rest of the 
potential. 

These models provide a very well-defined form for w{a) 
that depends only on the present-day values of w and Q^. 
While we have provided a variety of arguments both for 
and against such models, it is obvious that these issues 
will ultimately be settled by observational data, rather 
than by the speculations of theorists like us. 
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